Aggregation and fragmentation dynamics of inertial particles in chaotic flows 
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Inertial particles advected in chaotic flows often accumulate in strange attractors. While moving 
in these fractal sets they usually approach each other and collide. Here we consider inertial particles 
aggregating upon collision. The new particles formed in this process are larger and follow the 
equation of motion with a new parameter. These particles can in turn fragment when they reach 
a certain size or shear forces become sufficiently large. The resulting system consists of a large set 
of coexisting dynamical systems with a varying number of particles. We find that the combination 
of aggregation and fragmentation leads to an asymptotic steady state. The asymptotic particle size 
distribution depends on the mechanism of fragmentation. The size distributions resulting from this 
model are consistent with those found in rain drop statistics and in stirring tank experiments. 

PACS numbers: 05.45.-a, 47.52. +j, 47.53. +n 



There is an increasing recent interest in the advcction 
of inertial particles in fluid flows [1] . This comes in part 
from the fact that the dynamics of these particles are dis- 
sipativc. which leads in most flows to a preferential accu- 
mulation on chaotic, fractal attractors. Previous studies 
concentrated mainly on noninteracting particles, in spite 
of the fact that accumulation leads unavoidable to strong 
mutual interactions of different kinds. 

Here we consider the interaction of these particles in 
form of aggregation and fragmentation. When two par- 
ticles come sufficiently close, they aggregate to a larger 
one, observing mass and momentum conservation. If the 
size exceeds a certain threshold value, which depends on a 
property, the stickiness 7, of the particles, or other con- 
ditions arc fulfilled, they break up into smaller pieces. 
This is the basic mechanism underlying such processes 
in nature like raindrop formation in clouds [2] or the sed- 
imentation of marine aggregates [3] in the ocean. We 
demonstrate that for the study of such processes a parti- 
cle based approach may be a useful addition to the usual 
population-balance equation approach [4]. The latter is 
based on the assumption of well-mixed particles while our 
approach takes the incomplete mixing of inertial particles 
in fluids explicitly into account. 

Although concepts of dynamical systems theory can 
usefully be applied, we show that the entire dynamics 
is much more complex than that of any usual dynam- 
ical system. The dynamics of particles of any size are 
governed by the same type of equations of motion, but 
with different parameters since new particles will have 
new radii. Even if one considers a finite number, n, of 
possible sizes (size classes), there are n equations of mo- 
tion with different size-dependent parameters. We thus 
have a union of n dynamical systems and, moreover, the 
number of particles in each size class is changing in time. 
It is useful to interpret the attractors of the different size 
classes (of the non-interacting problem) as the skeleton of 
the full dynamics. Aggregation/fragmentation generates 



transitions from one attractor to another one. It is this 
permanent wandering among different attractors which 
characterizes the new dynamics. 

We show that the combination of aggrega- 
tion/fragmentation, superimposed on chaotic inertial 
advection dynamics, leads to a convergence to an 
asymptotic steady state, and this steady state is unique 
for the cases studied here. We find that the dynamics 
and the steady state depend on the fragmentation rule. 
For fragmentation due to shear we present a simple 
scaling relationship for the asymptotic average size of 
the particles. Furthermore, the shape of the asymptotic 
size distribution can be represented in a scaled form 
independent of the stickiness 7. 

For simplicity we consider spherical aerosoles, i.e. par- 
ticles much denser than the ambient fluid and assume 
that the difference between their velocity r and the fluid 
velocity u = u(r(i), t) at the same position is sufficiently 
small so that the drag force is proportional to this dif- 
ference (Stokes drag). The dimensionless form of the 
governing equation for the path Y{t) of such aerosols sub- 
jected to drag and gravity, reads as [5]: 



^(u 



(1) 



where n is a unit vector pointing upwards in the verti- 
cal direction. Throughout this paper we consider the 
vertical direction along the axis y. The inertia pa- 
rameter A (larger values for smaller particle size) can 
be written in terms of the densities pp and p/ of the 
aerosol and of the fluid, respectively, the radius a of 
the aerosols, the fluids kinematic viscosity v, and the 
characteristic length L and velocity U of the flow. It is 
A — R/ St, where i? = /9//pp-Clis the density ratio and 
St = {2a?U)/{9vL) is the so-called Stokes number of the 
aerosol [8]. W = 2a? ppg / {9v p fU) is the dimensionless 
settling velocity in a medium at rest. 

Every particle produces perturbations in the flow that 
decay inversely proportional to the distance from the par- 
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tide [6]. Here we assume a dilute regime, where the local 
concentration of particles is low enough, so that particle- 
particle interaction can be neglected [7]. 

During aggregation and fragmentation the radius of 
particles changes and so do the parameters A and W. 
The smallest (primary) particles considered in this model 
have dimensionless radius oi = 5/30^/"^ x 10~^, mass 
mi = /5p4/37raf , inertia parameter Ai = 7 and settling 
velocity Wi ~ QA/Ai. All larger particles are assumed 
to consist of an integer number of these primary parti- 
cles. 30 different size classes are considered. A particle 
that consists of a (a = 1, 30 is called size class index) 
primary particles has a radius Ga = a^^^ai, an inertia 
parameter A^ = (01/0^)^^1 = a~^^^Ai and a settling 
velocity = a'^^'^Wi. The largest particle therefore 
has a radius 039 = 5 x 10~^. 

Aggregation takes place upon collision, i.e. if two parti- 
cles, say of radius ai and aj, come closer than a threshold. 
Mass conservation requires the radius of the new particle 
to be a^g„, = a^+aj. For the size class index this implies 

a linear rule: anew = cti + oij which determines the new 
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inertia parameter via Aa^^^^ = anew ^i- The velocity of 
the new particle follows from momentum conservation. 

Fragmentation: we apply one of the following rules, (i) 
Size-limiting fragmentation: If a particle becomes larger 
than the maximum radius 030, it is broken up into two 
smaller fragments whose radii are chosen randomly, with 
a uniform distribution between oi and half the original 
radius. If any fragment is larger than 030 this process 
is repeated, until no fragment exceeds 030. (ii) Shear 
fragmentation takes place if the velocity gradient is too 
large. More specifically, the velocity gradient is evaluated 
across each particle both in the horizontal and in the ver- 
tical direction. If the maximum in any direction exceeds a 
threshold value, the particle is broken up into two smaller 
parts in the same way as for size-limiting fragmentation. 
While size-limiting fragmentation is dominant for rain- 
drops [2], shear fragmentation determines the break-up 
of marine aggregates [9] . 

Since for marine aggregates the threshold gradient be- 
comes smaller for larger particles [9], wc write 
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(grad(u))t,, = ^ai/a ja" 



(2) 



Coefficient 7 represents the 'stickiness' of the particles. 
Whatever rule is taken, the result is the reversed process 
of aggregation: two new particles are formed from an 
old one with the size class indices: ai^new + aj^new — 
aoid- The centers of the new particles are placed along a 
line segment in a random direction so that their distance 
equals the sum of their radii. Momentum is conserved. 
For simplicity we assume that the new particles have the 
same velocity as the old one. Shear fragmentation is 
applied together with size-limiting fragmentation to keep 
the maximum number of occurring size classes at 30. 

At the instant of both the aggregation and the frag- 
mentation process there is a sudden change in the dy- 



namics: the number of particles jumps in 3 among the 30 
available dynamical systems defined by the size classes. 

For convenience, we treat the case where the fluid flow 
is two-dimensional, therefore the phase space of the ad- 
vection dynamics is 4-dimensional. We use the convec- 
tion model of [10] with dimensionless velocity field 



u{x,y,t) ~ [1 + k sm{iLjt)] 



sin(27r.T) cos(27ry) 
— cos(27ra:) sin(27r?;) 



, (3) 



where k ~ 2.72 is the amplitude and uj = tt is the fre- 
quency of the periodic forcing. The fiuid flow itself is 
laminar, but the dynamics of the inertial particles can 
be chaotic. Because of the spatial periodicity of the flow 
and the resulting spatial periodicity of the attractors the 
total particle mass, M, in each 1x1 unit cell remains 
the same over time. The dynamics can therefore be re- 
stricted to one cell. The characteristic size and velocity 
of the flow are therefore L — 1, U = 1, respectively. 
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FIG. 1: Poincare section of the attractors of Eq. (1) projected 
onto the plane of the flow for inertia parameters (a) A = 7 
(size class 1), (b) A = 2.778 (size class 16), (c) A = 2.253 (size 
class 30). The positive Lyapunov exponents are (a) Ai = 
0.108, (b) Ai = 0.061, (c) Ai = 0.119, A2 = 0.014. The 
settling velocity is W — OA/ A. 

In the numerical realization of the problem the par- 
ticles are advected without any interaction over a time 
interval St = T/20 at the end of which first aggregation 
and then fragmentation take place, instantly. This is re- 
peated after every time step St. To carry out the aggrega- 
tion process, the distance between particles is calculated 
and all particles within a distance less than the sum of 
their radii aggregate. 

As initial condition we take 10^ particles in the smallest 
size class and no particles in other size classes. Further- 
more particles are uniformly distributed over the entire 
configuration space with velocities matching that of the 
fiuid. This choice fixes the total mass of the system to 
be M lO^mi. 

Before presenting the results obtained for the full dy- 
namics, it is instructive to see the attractors of the non- 
interacting problem. Fig. 1 presents the attractors for 
the smallest, an intermediate and the largest size classes. 
The extension of the attractor seems to grow almost 
monotonically with the size class index, except for a few 
intermediate size classes (a = 9. ..14), where the attractor 
size decreases or the attractor becomes periodic. 
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FIG. 3; (a) Average size class index {a{t)} vs. time for size- 
limiting fragmentation (upper curve) and shear fragmentation 
with 7 = 9 (lower curve), (b) Asymptotic average size class 
index for the same initial particle distribution as a function 
of the stickiness parameter 7. Squares; numerical results, 
continuous line: fit 000(7) by a linear function of 7'^, based 
on the data for 7 < 10 



FIG. 2: Particle numbers vs. time, and space distributions 
for size-limiting fragmentation, (a) Total number N{t) of par- 
ticles (bold - left axis) and the number of particles Na{t) in 
size class a (gray) for a — 1 (left axis), a = 16 (right axis) 
and a = 30 (right axis). Distribution of all particles in con- 
figuration space at time (b) t — 5, {c) t = 100. 



In order to understand the full dynamics, we include 
first the simplest fragmentation process, the size-limiting 
fragmentation. Fig. 2(a) shows the time dependence of 
the number Na{t) of particles in a few size classes. The 
particles leave the initial size class very quickly. After 20 
time units nearly all other size classes are considerably 
occupied. In fact, the population in size class 16 reaches 
a maximum here, but decreases again later on. It is the 
occupation of the largest size classes which continuously 
increases and then saturates. The total number N{t) 
of particles (bold line) rapidly decreases first, but satu- 
rates later on. The spatial distribution of particles (Fig. 
2(b/c)) shows that they move initially among the more 
localized attractors characteristic of small size class in- 
dices. Later, the distribution becomes more extended in 
configuration space when size classes with extended at- 
tractors become well occupied, although the total number 
of particles is much less than in the initial phase. While 
the full dynamics is dominated by transients in-between 
attractors, the shape of the backbone attractors is clearly 
recognizable in the plots. 

To follow the convergence towards an asymptotic state, 
we found useful to consider the average size class index 
{a{t)) = 'E^'Lia^Na^{t)/N{t). Fig. 3(a) shows the time 
dependence of this index for both types of fragmentation. 
It illustrates the convergence to an asymptotic steady 
state for both fragmentation rules. Initially, aggregation 
leads to a fast increase in the average particle size class 
for both fragmentation rules. Then fragmentation sets 
in and a balance between aggregation and fragmentation 
is reached, with a different asymptotic average particle 



size Q!oo — lim {a{t)) for the two rules. For size-limiting 

t — >oo 

fragmentation the value of a^^ is almost constant over 
time, while for shear fragmentation a^o oscillates with 
the period T of the flow. This is caused by the periodic 
change in the fluid flow and the corresponding change in 
the shear forces. 

For size-limiting fragmentation, ctoo is, in a broad 
range, independent of M. For shear fragmentation with 
M < 3 X lO^mi, aao{M) increases approximately linearly 
with M, while for higher values a saturation of aoo{M) 
sets in, which is due to size-limiting fragmentation. 

By considering other initial conditions than mentioned 
above, while keeping the total mass M fixed, the asymp- 
totic state is found for both rules to be independent of 
the chosen initial condition, but for shear fragmentation 
the asymptotic state does depend on the value of the 
stickiness 7. 
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FIG. 4: Histogram of the particle size distribution, (a) Size 
class percentage Na/N in the steady state versus the dimen- 
sionless radius a, for size-limiting and shear fragmentation 
with 7 = 9. Results for fragmentation into three parts is 
also shown (open markers), (b) Normalized number of parti- 
cles versus the relative radius a/ (a) for different values of 7 
(binary splitting). 

To illustrate this dependence of the steady state on the 
stickiness 7, Fig. 3(b) shows how aoo changes with the 
stickiness parameter, at a fixed M. A drastic increase 
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of a(7) can be observed in the interval 4 < 7 < 10. It 
is clear that aoo increases with 7, because particles be- 
come more resistant to shear. A quantitative estimate of 
the shape of this (7) curve can be derived by assum- 
ing that the threshold velocity gradient is approximately 
constant for the size class index aoo- From (2) it then 
follows that aoo depends linearly on 7'^. This simple de- 
pendence is expected to hold for relatively small values 
of 7 and aoo, where shear fragmentation dominates. It 
can be seen that for higher values of 7, when size-limiting 
fragmentation becomes important the aoo (7) curve devi- 
ates from this estimate and converges towards a limiting 
value a^™^ (Fig. 3(b)). 

In addition to the average quantities it is natural to in- 
vestigate the occupation of the different size classes in the 
steady state. Fig. 4(a) shows the steady state histograms 
vs. the dimensionless radius. For size-limiting fragmen- 
tation the distribution shows one broad peak around 
smaller size classes and a second, smaller peak at large 
size classes with a sharp drop-off towards zero beyond 
the maximum size. This behavior, with two maxima and 
a sudden drop after the second peak, is similar to that 
of observed cloud drop spectra [2] . For shear fragmenta- 
tion the steady state distribution also shows two peaks, 
but much closer together, with a long tail in the parti- 
cle distribution towards larger sizes that goes smoothly 
towards zero. For 7 < 5 this distribution is not fully 
developed and only shows one peak. In the intermedi- 
ate 7 range, where the distribution is fully developed, 
but size-limiting fragmentation is not important, a scal- 
ing form = (w) iound {{a) represents the 
average radius), independently of 7. All distributions in 
this range collapse then onto a single master curve as 
shown in (Fig. 4(b)). This behavior, along with the long 
tail in the distribution towards the right hand side, is 
typically observed in shear-fragmentation experiments in 
stirring tanks [11]. 

We note that our findings are robust with respect to 
the number of new particles formed by fragmentation. 
For instance, in Fig. 4(a) we see that the distributions 
of particles for ternary fragmentation are similar to the 
ones for binary splitting and only show a slight shift to- 
wards smaller size classes. The same result is found in 
population balance equation models, e.g. [11]. 

Finally wc mention that in spite of the different steady 
states, the size distribution in the initial phase is simi- 
lar in the different cases. After short times, wc find a 
roughly exponential decay. In this early phase, fragmen- 
tation is yet inactive, and the process is dominated by 
aggregation. This decay in the short time distribution 
can be found for all initial conditions. 

In conclusion, we illustrated that an individual model- 
ing of particles is able to reflect typical properties of ag- 



gregation/fragmentation processes. We found the devel- 
opment of a balance between aggregation and fragmen- 
tation, and a steady state. The steady-state particle size 
distributions found here correspond to those observed in 
rain drops (size- limiting fragmentation) and stirring tank 
experiments (shear fragmentation). For shear fragmen- 
tation the size distributions are found to follow a scaled 
form. In addition the approach shown here can reflect 
spatial inhomogeneity and take actual particle dynam- 
ics into account, and could possibly allow for a much 
more detailed description of particle interaction. It is 
thus more adequate than the usual stochastic, mean field 
like approach which relies on the assumption that the 
particles are well mixed [2]. The presence of chaotic at- 
tractors can ensure a partial mixing only and hence the 
assumption is not valid. 

An interesting open problem is to extend our study to 
3D flows. 
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